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Abstract 

We calculate the geometric phase associated with the time evolution of the wave function of 
a Bose-Einstein condensate system in a double-well trap by using a model for tunneling between 
the wells. For a cyclic evolution, this phase is shown to be half the solid angle subtended by 
the evolution of a unit vector whose z component and azimuthal angle are given, respectively, 
by the population difference and phase difference between the two condensates. For a non-cyclic 
evolution, an additional phase term arises. We show that the geometric phase can also be obtained 
by mapping the tunneling equations on to the equations of a space curve. The importance of a 
geometric phase in the context of some recent experiments is pointed out. 
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I. INTRODUCTION 



Bose-Einstein condensation (BEC) in a dilute gas of trapped ultracold alkali atoms has 
been observed by several experimental groups [1] . This gives rise to the possibility of under- 



standing the nature of the condensate wave function, and in particular, its phase [2j. It is 
believed that the BEC phase transition occurs due the breaking of a global gauge symmetry 
of the Hamiltonian. Theoretically, a BEC may be modeled by writing down the interact- 
ing many-body Hamiltonian in terms of boson creation and annihilation operators ^ and 
ty op . The order parameter is postulated to be the condensate wave function ip—< ty op >= 
pe , where p=\ip\ 2 is the condensate density and 9 is the phase of the wave function. The 
Hamiltonian is gauge invariant, but the order parameter breaks this symmetry. Using the 
dynamical equation for ^> op found from the Hamiltonian operator, the time evolution of the 
condensate wave function if> can be shown to satisfy the following Gross-Pitaevskii equation 
(CPE) S: 

%K ^t = ^ ^ + ^ t(x) + 9o ^ ^ (L1) 

where V ex t is the external potential and go = 4:nh 2 a/m, a and m being the atomic scattering 
length and mass, respectively. Although this equation has an underlying quantum nature, 
the condensate has a macroscopic extent, suggesting the observation of quantum effects on 
a macroscopic scale. 

In a striking experiment Q , Andrews et al have shown the existence of the macroscopic 
quantum phase difference between two bose condensates : They designed a double-well trap 
by using a laser sheet to create a high barrier within a trapped condensate. On switching 
off this barrier, the two condensates overlapped to produce an interference pattern, showing 
phase coherence. More interestingly, by lowering the laser sheet intensity, the barrier gets 
lowered, making it possible for the condensates to tunnel through the barrier. Thus this 
double-well trap is analogous to a superconductor Josephson junction p, and is referred to 
as the Boson Josephson junction (BJJ). In an interesting paper, Smerzi et al [6j have set up 
the tunneling equations for the BJJ in a model. These are two coupled nonlinear ordinary 
differential equations for the condensate wave functions in the two wells. They have studied 
the time evolution of the inter-well population difference and phase difference in this model, 
and predicted a novel 'self-trapping' effect, i.e., the oscillation of the population difference 
around a non-zero value, for certain initial conditions and parameters. 
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The tunneling dynamics motivates the following question: Is there an underlying geomet- 
ric phase associated with the time evolution of the condensate wave function in a double- well 

trap ? As is well known by now, the concept of a geometric phase has been studied in vari- 

n 

ous contexts, after it was introduced by Berry |2j in quantum mechanics. It had also been 
considered much earlier by Pancharatnam [^| in the context of classical optics. Geometric 
phase and its various applications have been studied intensively for over a decade now 
Such a (non-integrable) phase arises when the time evolution of a system is such that the 
value of a variable in a given state of the system depends on the path along which the state 
has been reached. In this paper, we calculate the geometric phase associated with the time 
evolution of the BJJ wave function, for both cyclic and non-cyclic evolutions. 

The plan of the paper is as follows: In Sec. II, we first review Smerzi et al's [f| derivation 
of the tunneling equations. Keeping in mind that the geometric phase is gauge independent, 
we use a certain gauge transformation to reduce these equations to a more convenient form, 
n Sec. Ill, we briefly outline the kinematic approach formulated by Mukunda and Simon 
lfll | to define the geometric phase as applied to a two-level system. We then solve the BJJ 
tunneling equations, which are nonlinear differential equations, numerically by choosing 
some parameter values as an example. Using these solutions, we find the geometric phase 
explicitly, for both cyclic and non-cyclic evolution of the system. For a cyclic evolution, 
the phase difference and the population difference between the condensates in the two wells 
return to their original values. For this case, the corresponding geometric phase is half the 
solid angle generated by a unit vector whose z component and azimuthal angle are given, 
respectively, by the population difference and the phase difference. For non-cyclic evolution, 
an additional phase term is obtained. In Sec. IV, we show that this geometric phase (for 
both types of evolution) can also be obtained by first mapping the tunneling equations to 
the equation for a unit vector r and then identifying it with the tangent T of a space curve. 
The space curve is described using the so-called natural frame equations, which possess an 
underlying natural gauge freedom. The unit triad of vectors can be written down using 
the form of the condensate wave functions in the two wells. The concept of Fermi- Walker 
parallel transport is then used to identify the geometric phase. In Sec.V, we employ the 
usual Frenet frame to obtain explicit expressions for the curvature and torsion of the space 
curve that gets associated with the BJJ evolution. Section VI contains a summary and 
discussion. 
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II. THE BJJ TUNNELING EQUATIONS 

We begin by briefly describing the model used by Smerzi et al6] to study the tunneling 
of the condensate between two wells. Let the total number of atoms in the double-well trap 
be N. Let Ni and N 2 denote the number of atoms in each well, such that Ni + N 2 = N. To 
study the tunneling, the solution for the GPE (Eq. (1.1)) is assumed to be of the form 

V> = Vi +V> 2 (t)$ 2 (x). (2.1) 

Here $i,$2 are the ground state solutions for the isolated wells with Ni = N 2 = (N/2). 
Using Eq. (2.1) in Eq. (1.1), one gets 

ih^ = (£? + C/xiVOVr - Vijj 2 (2.2a) 

ih^ = (£ 2 ° + U 2 N 2 )ip 2 - Vfa. (2.26) 

The quantities E® 2 ,Ui t2 and V are constants obtained by taking the external field V ext (x) 
in Eq. (1.1) to be independent of time: 



K2 = J [(ft 2 /2m)| V $i, 2 | 2 + $?, 2 V ext ) d 3 x 
Ui )2 = 9o J ®i 2 d 3 x, 



y = -J i(h 2 /2m)( V $i V $2) + $i$2 Vest] d 3 x 

Further, Ni = \ipi\ 2 , for i — 1,2. 

In this paper, we obtain the geometric phase associated with the evolution in Eqs. (2.2). 
Recognizing that this phase is a gauge independent quantity, we find it convenient to apply 
the following gauge transformation to the wave function 




JfeSvW* I a I . (2.3) 



On substituting Eq. (2.3) in Eq. (2.2) and setting hrj(t) = — ((J5° + E$) + N{U x \a\ 2 + 
^ r 2|6| 2 ))/2, the BJJ tunneling equations take on the form 

■\L,. I , (2.1) 





with 



hw = - El + UiN x - U 2 N 2 ) (2.5) 



In Eq. (2.4), M Wo denotes the (time-dependent) matrix on the right hand side. The 
advantage of the gauge transformation is that this matrix is traceless. Note also that since 
Ni = N\a\ 2 and N 2 = N\b\ 2 appear (in uq) in this matrix, Eqs. (2.4) are two coupled 
nonlinear differential equations. From Eq. (2.3), we see that the normalization condition 
l^i 1 2 + I V*2 1 2 — N implies |a| 2 + |o| 2 = l. Thus without loss of generality, we write 

a = cos(a/2) e Wl ; b = sin(a/2) e id2 (2.6) 

Let us denote the difference in the population density of the two traps by z and the difference 
in the phases of the two condensates by 0. From Eq. (2.6) we thus have, 

z = (N 1 -N 2 )/N=(\a\ 2 -\b\ 2 ) = cosa ; <j>=(0 2 -O 1 ). (2.7) 

By suitably combining Eq. (2.4) and its complex conjugate, and using Eq. (2.7), the 
nonlinear coupled equations for z and are found to be (on setting h — 1) 

_ = -Vy/1 - z 2 sin (2.8a) 

j = Az + V I cos + AE. (2.86) 
n at vl-2 2 

Here we have defined |6| 

AE = {El - El)/2 + N{U X - U 2 )/A ■ A = N{Ut + U 2 )/A (2.9) 

and the time has been reparametrized as t — * 2t. It is interesting to note that the above 
equations can also be written as Hamilton's equations, by treating z and as the canonically 
conjugate variables. The classical Hamiltonian is easily verified to be 

H cl = A VVl-z 2 cos + AEz. (2.10) 

2 



This describes a non-rigid pendulum with a length proportional to yl — z 2 , which decreases 
with the "momentum" z. We write down the expression for uq appearing in Eq. (2.5) in 
the form 

u = AE + Az, (2.11) 
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by using Eq. (2.9). Finally, setting z = cos a in Eqs. (2.8a) and (2.8b), we obtain 




(2.12a) 



— = A cos a + V cot a cos <fi + AE 



dt 



(2.126) 



Equations (2.8), or equivalently, Eqs. (2.12), represent the tunneling equations. For 
convenience, V can be absorbed in time t and all energies can be measured in units of V. 
Note that AE is the asymmetry between the two wells, as seen from Eq. (2.9). We consider 
two special limits: 

(1) Interacting Bose system in a symmetric trap: A^O, AE = 0. 



From Eq. (2.9), £° = E° and U x = U 2 = U, and A = UN T is nonzero. Thus Eq. (2.10) 



In Figs. 1, 2 and 3, we have obtained the (z,<j>) phase portraits for this case, by using 
Eq. (2.13). Let A be replaced by the dimensionless quantity (A/V). For A < 1 there exist 
periodic oscillations around the zero-state (0,0) and the non- trapped ir state (0, 7r). There 
are no rotational states, (see Fig. 1 for A = 0.5). For A > 1 two new trapped 7r-states 
appear at (z*, ir) and (—z*, ir) with z* = VA 2 — 1/A. The traped ^-states are clearly visible 
in Fig. 2 (A = 1.3). As A increases, z* — > 1. For A > 2, rotational states also appear as 
seen in Fig. 3 (A = 5). 

(2) Non-interacting Bose system : A = 0. 
For an ideal bose gas, the interactions U\ and Ui are negligible, and Eq. (2.9) yields A = 0. 
In this limit the kinetic energy term of the Hamiltonian H c \ (2.10) vanishes and hence the 
non-rigid pendulum analogy is not valid anymore. (However, we remark that the tunneling 
Hamiltonian coincides with that of a two-component BEC in the rotating frame approxima- 
tion. This will be discussed in the last section.) We get 



In Figs. 4, 5 and 6, we have obtained the phase space portraits for this case. Again, as in 
case (1), the energy is measured in units of V. When AE = 0, we have a symmetric trap. 
There are oscillations around the (0, 0)-state and the non-trapped (0, 7r)-state. There are no 



becomes 




(2.13) 



H d = -VVl - z 2 cos + AEz 



(2.14) 
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rotational orbits, as seen in Fig. 4. With increase in AE, rotational orbits appear, which 
explore the full range of 0. This is shown in Figs. 5 and 6. The oscillations are now around 
(z = —z*, = 0) and (z = +z*, = ±7r). That is, oscillations around = shifts towards 
z — — 1 which is energetically more stable, while the = ir fixed point moves towards z = 1 
which is a local energy maximum. Note that the increase in AE causes a non-trapped ir 
state to become a trapped 7r state with oscillations around a non-zero population difference 
z*. The new fixed point is given by z* = \J '((AE /V) 2 / '(1 + (AE/V) 2 )) which tends to 1 as 
AE/V goes to infinity. 

We now proceed to show how the geometric phase associated with the BJJ dynamics can 
be computed. 

III. GEOMETRIC PHASE USING THE KINEMATIC APPROACH 

In this section we derive the expression for the geometric phase for the BJJ evolution 
using the kinematic approach developed by Mukunda and Simon 10]. We first briefly outline 
the basic ideas. One starts with a complex Hilbert space and considers a subset Ao made up 
of unit complex vectors denoted by ip(t). Let Cq be a smooth one parameter curve, consisting 
of a family of vectors ip(t), with < t < T . Since ip(t) are unit vectors for all values of t, 

ife(V>(*),#)) = o, 

giving 

m)M)) = ^m),m) (3.1) 

where the dot denotes the derivative with respect to t. Under a gauge transformation defined 
by a real function a(t), C -> C and if)(t) -> ^'(t) = e ia ^i)(t), t £ (0,T), one finds, 

1mW(t)rf(t)) = lm(m,m) + (3-2) 
On the other hand, we have 

arg(i;'(0)^'(T)) - arg(ip(0),ip(T)) = dt 
Integrating Eq. (3.2) and using the above equation, we see that a gauge invariant function, 

arg(ij'(0),ij'(T))-lm f dt (^'(t)^(t)) (3.3) 

Jo 
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can be constructed. Next one introduces the space TZo of unit rays, which is the quotient 
of Ao with respect to the U(l) action ip — ► e ia vp. It can then be shown that the curves 
C and C' (obtained by a gauge transformation) in Ao both project to the same curve c 
in the ray space 1Z . Thus the gauge invariance of the functional (3.3) implies that it is a 
functional only of the image Cq of Co- This functional can be shown to be reparametrisation 
invariant as the integrand is linear in ip'(t). One thus defines the geometric phase 9 (co) as 
the following gauge and reparametrisation invariant functional, 

9 [c o ] = an?(V/(0), V'(T)) - Im f dt ^\t)^\t)) (3.4) 

Jo 

Identifying the first term in Eq. (3.4) as the total phase 

<P P = arg(^(0)^'(T)) (3.4a) 
and the second term in Eq. (3.4) as the dynamical phase 

</> d = lmf dt (^(t),^'(t)), (3.46) 
Jo 

we get 

4> g [co\ = 4>p[C } - <j)d[C Q ], (3.5) 

Now for a given c G TZq the gauge freedom in the choice of C G Ao can be used to 
express 9 [co] in different but equivalent ways. The Pancharatnam connection takes ip(0) 
and ip(T) to be in phase, giving <p p = 0, leading to <p g = —<pd- Another prescription takes a 
lift of Co so that the integrand of (pd in (3.4b) vanishes. This leads to <p g = <p p . 

Let us now understand what these two prescriptions imply for the BJJ evolution equa- 
tions. Using equation (2.6) the family of unit vectors is given by, 



</>'= 1 




| = e">>«> ( 




^1 





cos(a(*)/2) 
sin(a(t)/2)e^W 



= e i9l(t V, (3.6) 



where = (9 2 — 9i). Clearly, the angle 6\ in Eq. (3.6) represents the gauge freedom in the 
wave function ip'(t). The prescriptions essentially fix the gauge 9±. Using Eq. (3.6) in Eq. 
(3.4a), a short calculation leads to the following total phase: 

P = arg(^(0)^'(T)) = (9 1 (T) - 9,(0)) + A. (3.7) 

Here, 



A = tan" 1 



sin(a(0) /2) sin(a(T) /2) sin(0(T) - 0(0)) 



(3.8) 



_cos(a(0)/2) cos(a(T)/2) + sin(a(0)/2) sin(a(T)/2) cos(0(T) - 0(0)). 

The integrand of the dynamical phase 0^ can be calculated using Eq. (3.6) in Eq. (3.4b) to 
give 

Im(^', ( ^-) = 6>i + sin 2 (a/2)0. (3.9) 
at 

Thus 

pT 

(3.9a) 



(3.96) 



<f> d = (9 1 (T)-e 1 (0))+ / sin 2 (a/2)0 dt. 

Jo 

From Eqs. (3.7) and (3.9a), we get the geometric phase to be 

<p g = <p p - 4>d = - I sin 2 (o;/2)0 dt + A. 
Jo 



Next, we also calculate the geometric phase expression for two different gauge prescriptions, 
for the sake of completeness. 

The Pancharatnam Prescription. 

Here P = 0, thus Eq. (3.7) yields, 



e 1 {T)-e 1 (o) = -A. 



(3.10) 



Given the end-points a(0), a(T), 0(0), 0(T) on the curve C one can thus determine A from 
Eq. (3.8), and hence the gauge $i{T) - 0i(O). From Eq. (3.9), 



<p d = Jlm(ij},i}})dt = e 1 (T)-e 1 (0) + J sin 2 (a/2)0dt. 
Substituting Eq. (3.10) in Eq. (3.11), 

4> d = -A + / sin 2 (a/2) 0dt. 
Jo 



(3.11) 



Since 



i d , we get 



f sin 2 (a/2)0 dt + A, 
'o 



(3.12) 



(3.13) 



which agrees with Eq. (3.9b). 
The Horizontal Lift. 

Here the integrand in the 0^ definition is zero, giving 



0i = - sin 2 a/2 0, 
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so that the gauge is determined in this case to be 

0i (T) - 9 1 (0) = - I sin 2 (a/2) di. (3.14) 

Substituting this in Eq. (3.7), 

<f>p=- f (sin 2 a/2) dt + A. (3.15) 

Since 0^ = here, we get an expression for 9 which is the same as Eq. (3.13). Thus the 
two prescriptions give the same geometric phase, which holds for both cyclic and noncyclic 
evolutions. 

For a cyclic evolution, it is clear from (3.8) that A = 0. Hence the geometric phase is 
just minus half the solid angle Q subtended by the closed curve generated on a sphere by 
the tip of a unit vector r: 

r = (sin a cos 0, sin a sin 0, cos a) (3.16) 

Here, a and denote the polar and azimuthal angles of r. 
An Example: 

As an example we consider an interacting Bose system withA = 5, in a symmetric trap, 
i.e., AE = 0. As already mentioned, the phase space portrait for this case is found from 
Eq. (2.13) and is given in Fig. 3. The value of A selected is quite generic because further 
increase in A does not change the character of the phase- space portrait much, apart from 
moving the 7r-state towards z — 1. Since cos a = z, the geometric phase <p g given in Eq. 
(3.9b) can be re-expressed as 

<t> g = \ I (z-l)j>dt + A. (3.17) 
^ Jo 

While comparing the geometric phase from the above equation to that in Eq. (3.4), it is 
necessary to keep in mind that in the Hamiltonian formalism the time has been scaled and 
so one has to use the appropriate value of time in Eq. (3.17). We solve Eqs. (2.8a) and 
(2.8b) numerically for (z(t) , (f>(t)) , for a given initial condition (z(0), 0(0)) at time t — 0. 
Using these solutions, we find the solutions for the corresponding unit vectors r given in Eq. 
(3.16). These yield the the path on the unit sphere plotted in Fig. 7. Next, we substitute 
the solution for (z(t),<f>(t)) in Eq. (3.17), to find 4> g {t) numerically, for various times t, in 
the range < t < T, where T denotes the full period of the orbit concerned. Our plot for 
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the time dependence of <p g (t) over a period for a librational orbit is given in Fig. 8, while 
Fig. 9 gives that plot for a rotational orbit. 

We conclude this section with the following remark. There exists an interesting geomet- 
rical representation of a two level system in terms of the time evolution of a unit vector r. 
In the next section, we identify r with the tangent of a space curve, and provide a classical 
differential geometric approach to derive the geometric phase <p g associated with the BJJ 
evolution. 



IV. GEOMETRIC PHASE USING SPACE CURVE APPROACH 



In this section, we derive the geometric phase associated with the Bose condensate tun- 
neling dynamics by providing a geometric visualisation of this two level system. Firstly it 



is possible to show 



111 ] that the tunneling equations (2.4) for the two-level wave function, 



( 













cos(o;/2) 
sin(a/2)e^ 



with = 6 2 — Q\ can be mapped to the following vector evolution equation. 

dr/dt = u) x r. 



Here, in cartesian coordinates , 



(4.1) 



(4.2) 



u = (-2V, 0, 2u n 



(4.3) 



r = (a*b + ab*, i{ab* - a*b), \a\ 2 - \b\ 2 ). (4.4) 

Using the definitions of a and b given in Eq. (2.6), Eq. (4.4) is readily seen to be identical 
to the unit vector r in Eq. (3.16). 

has shown that corresponding to a wave function of the form (4.1), in 
addition to r two more unit vectors P', Q' can be defined, such that the set (r, P', Q') forms 
a unit orthogonal right-handed triad. This is achieved by defining a complex vector Z' as 
follows : 

Z' = P' + iQ' = ((a 2 - b 2 ), i(a 2 + 6 2 ), -2ab). (4.5) 
On using Eq. (2.6) in Eq. (4.5), we get, 
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71 = e 2iei (cos 2 (a/2) - sm 2 (a/2)e 2 ^, i(cos 2 (a/2) + sin 2 (a/2)e 210 ), - sin(a)e^) (4.6) 

It can be easily verified that, 

|r| = |P'| = |Q'| = 1 ; r-P' = r-Q' = P'-Q' = (4.7) 

Clearly now as r evolves with time, so does the (P', Q') plane. The total phase T p 
accumulated by Z'(£) in time T is given by, 

F p = arg(Z'(0)* ■ Z'(T)) (4.8) 

Substituting for 7J from Eq. (4.6) into Eq. (4.8), after some algebra we obtain, 

r p = 2[(0 1 (T)-0 1 (O)) + A] (4.9) 

where A is identical to the expression Eq. (3.8) obtained in the kinematic approach in Sec. 
III. Thus from Eqs. (3.7) and (4.9) we get 

T p = 20 p (4.10) 

(f) p being the total phase of vp in the kinematic approach. 

First we find the total phase rotation 7 p associated with the rotation of P' or (Q') as 
follows. It is defined by 

cos 7p = P'(T) • P'(0) = Q'(T) • Q'(0). 

Further, it is easy to see geometrically that P'(T) • Q'(0) = — Q'(T) ■ P'(0) = sin7 p . Substi- 
tuting Z' = P' + iQ' in Eq. (4.8) and using the above relations, we can show that the total 
phase 

7p = -r p = -2[(^(T) - 6 1 (0)) + A], (4.11) 

where we have used Eq. (4.9). 

Next we wish to find the dynamical phase 7^ associated with (P', Q') rotation, which is 
induced by the specific dynamical equations of the frame (r,P', Q'). This is a little more 
involved, and we proceed as follows. 

From Eq. (4.6), we have 

Z' = e 2 ' Wl Z (4.12) 
12 



This immediately leads to 

P' + iQ' = e 2l9l (P + iQ), (4.13) 
Comparing this with Eq. (4.6) yields 

P = (cos 2 («/2) - sin 2 (a/2) cos 20, - sin 2 (a:/2) sin 20, - sin a cos 0) (4.14a) 

Q = (-sin 2 (a/2) sin 20, cos 2 {a/2) + sin 2 (a/2) cos 20, — sin a sin 0) (4.146) 

It can be easily verified that (r, P, Q) is also a right-handed triad. 

A short calculation using Eqs. (3.16) and (4.14) shows that we can write 

dr/dt = XP + YQ, (4.15) 

where 

da d(f) 
X = (— ) cos0 — (since— ) sm0 (4.16a) 
dt dt 

d(f) da 

Y = (sin a— ) cos0 + (— ) sin0 (4.166) 

(XL (JjL 

Obviously, there is a gauge freedom 29i in the choice of (P, Q). We immediately see this 
from Eq. (4.13): 

P' = Pcos/3-Qsin/3 (4.17a) 

Q' = P sin /3 + Q cos (3, (4.176) 

where 

13 = 20i (4.18) 

represents the gauge freedom. Using Eqs. (4.17), we solve for (P, Q) in terms of (P', Q'). 
Substituting them in Eq. (4.15) yields 

dv 

— = a l P' + a 2 q! (4.19) 

where 

ai = -7- cos(0 + 13)- (sin a^j-) sin(0 + (3) (4.20a) 

(JjL (JjL 

a 2 = ^- sin(0 + 13) + (sin a^-) cos(0 + (3) (4.206) 
dt dt 

Since (r, P', Q') is an orthonormal triad, Eq. (4.19) immediately implies, 

dP' 

— = -a 1 r + a 3 Q' (4.21) 
13 



at 

where a 3 is to be determined. In the space curve language, if r is identified with the 
tangent T, then a.\ and a 2 are the components of the curvature vector ^ along P' and Q' 
respectively. Further, equations (4.19), (4.21) and (4.22) describe the equations for a space 
curve in a "natural frame" (T, P', Q'). We remark that the Frenet frameQ corresponds to 
«2 = 0, P' is the normal n, Q' is the binormal b. Further, a 3 is the torsion r and a,\ is the 
curvature K. On setting a 2 = 0, we get from Eq. (4.20), the following "Frenet gauge" (3p: 

tan(/? F + 0) = sina^/A (4.23) 
at at 

Working with the natural frame, a short calculation using Eqs. (4.15) to (4.18) yields, 

a 3 = T • (t x T)/|T| 2 - 4- tan- x (— ). (4.24)(4.24) 

dt d\ 

Next using the cartesian representation of T = r given in Eq. (4.4), a lengthy but 
straightforward calculation leads to, 



T • (t x T)/|T| 2 = cosa^ + ^-tan^ 1 

dt dt 



sin a-jr 

dt_ 

da 
dt 



(4.25) 



Substituting Eq. (4.25) and (4.20) in Eq. (4.24) and using the formula tan 1 A — 
tan' 1 B = tan" 1 ((A - B)/(l + AB)), we obtain, 

dcj) d(<j> + 0) . 2 otd<f> d(3 

a? = cosa— = — 2 sin — — 4.26 

dt dt 2 dt dt v 1 

Note that the time derivative of the gauge freedom (3{t) appears in a 3 . 
We write Eqs. (4.19), (4.21) and (4.22) in a compact form, 

£= txTi f =«XP<; f={xQf. (4.27, 

Here $, is given by, 

| = « 3 T + aiQ' - a 2 P'. (4.28) 

Eqs. (4.27) show that the natural frame (T, P', Q') rotates with an angular velocity £, 
as it moves along the space curve. As is obvious, ax and c%2 are components of £ along the 
Q' and P' axes respectively and hence tilt the (P', Q') plane. On the other hand, a 3 merely 
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rotates^ plane around T. Thus in time T, the (P<,Q<) plane gets rotated bvan angle 
Id = J ot^dt. Such a frame is defined using Fermi- Walker parallel transport as (14J , 

= {( ai Q'-a 2 P')xA}\ 

dt 

Using the expression for as given in (4.26) we obtain the dynamical phase 7^ associated 
with (P', Q') plane to be 

ld = £ a 3 dt = -2 ^ T (sm 2 |) ^ dt - 2(9 1 (T) - 9 1 (0)) , (4.29) 

since (5 = 26 h from Eq. (4.18). 

Subtracting Eq. (4.29) from the expression for the total phase 7 P given in Eq. (4.11) we 
obtain the geometric phase 7 5 associated with (P', Q') rotation to be 

7g = 7 P -7, = 2[^ T (sin 2 |)^ dt — A] (4.30) 

Note that the term involving the gauge freedom f3 cancels out here too, as in the kinematic 
approach. Comparing Eq. (4.30) with Eq. (3.9b), we see that 

In other words, the geometric phase <p g associated with the wave function is minus half 
of that associated with the (P', Q') rotation. For a cyclic evolution, A = 0. Here, on 
computing j g , the geometric phase (fi g becomes just minus half the solid angle, as we saw 
in Sec. III. In summary, by mapping the evolution equation for the wavefunction to the 
dynamical equation for an orthonormal triad (T, P', Q') and identifying the triad to be a 
natural frame on a space curve, enables us to provide a purely geometrical visualisation of 
the geometric phase of a two level system. 

Our general result is valid for any two level system with the wavefunction (4.1), since 
we did not use the specific BJJ equations (4.2) and (4.3) in its derivation. By finding 
the solutions a(t) and <p{t) for the nonlinear equations (2.12) numerically for given initial 
conditions, 7 9 can be computed and is exactly -2<f) g , with <p g values as plotted in Figs. 8 and 
9. 



V. GEOMETRIC PARAMETERS ASSOCIATED WITH BJJ DYNAMICS 



In the last section, we discussed the mapping of the BJJ tunneling equations to a space 
curve which is described using equations for a "natural frame". This description involves 
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three geometrical parameters a% which are shown to depend on a gauge parameter (3 (see 
Eq. (4.20) and (4.24)). 

The usual description of a space curve is in terms of a Frenet frame [l^ . with the curvature 
K and torsion r as the geometric parameters. As explained in Sec. IV, working with the 
Frenet frame implies fixing (5 = /3p, defined in Eq. (4.19). In this section we work with the 
Frenet frame to determine the geometric parameters K and r of the space curve associated 
with the BJJ dynamics, in terms of the physical parameters V, AE and A and discuss certain 
special cases of interest. 

As mentioned in Sec. IV, in the Frenet frame, ai = K, a% — 0, = r, P = n and 
Q = b in Eqs. (4.15) to (4.18). In this frame, we have the usual Frenet-Serret equations 

— = Kn , — = -AT + rb ; — = -rn (5.1a) 

dt dt dt v 1 

Thus, 

dT . _ dn dh ... 

*- = e ' xT ; *=^ xn ; *=«- xb - (5 - 16) 
K2 = { it )2= sin2a( ^ )!+( ^' 2 (5 ' 2) 

•■ x , 9 dd> d i ,siaaf , 

r = T ■ (T x T)/A 2 = cosa^- + - tan"^— ^-), (5.3) 

dt 

where the cartesian representation (3.16) , T = r has been used. On using Eqs. (2.12) for 
^, ^ in Eqs. (5.2) and (5.3) respectively, we get 

K = 2{V 2 + K 2 ujI sin 2 a -V 2 sin 2 a cos 2 <p 

+2\^/io;oCOSQ;sinQ;cos0) 1 ^ 2 , (5.4) 

and 



Also, 



and 



t = cos a{hjjQ + V cot a cos (p) 

d _i .huo sina + V cos a cos d) _ N 

+ -i-tan ( , ). (5.5) 

ar 1/ sm <fi 

Equations (5.4) and (5.5) give the curvature and torsion of the space curve created by the 
BJJ dynamics Eq. (2.4), with parameters V and Jtloq. Here a and <fi are solutions of the 
integrable equations (2.12). 
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Since r is identified with T, we also have, for the BJJ system, 



dT 

— =u>xT, u> = (-2V,0,2a;o), (5.6) 

From Eq. (4.2) therefore K 2 can also be written as, 

K 2 = (^) 2 = M 2 - (« ■ T) 2 . (5.7) 

Using the definition of the matrix M Wo given in equation (2.4), a simple calculation shows 
that 2 < M wo >= uj • T, yielding 



4 < Ml o >= (a;) 2 . (5.8) 
Using Eq. (5.8) in Eq. (5.7), we get 

K = 2(< M 2 , > - < M W0 > 2 ) 1 / 2 . (5.9) 

Now from the first equation in Eq. (5.1) it is clear that the distance traveled by the tip 
of T on the unit sphere in time dt is ds = Kdt.This is the well known ^| Fubini-Study 
metric. Thus we see that the curvature K which determines the geometric quantity ds is 
given by the variance of the tunneling matrix M m (see Eq. (2.4)) for a two level system. 
As seen from the equation (5.3), the torsion integral J rdt measures the anholonomy of the 
frame, i.e. a path dependent geometric quantity given by the solid angle associated with a 
cyclic evolution of T. 

Recalling that the population density difference between the two traps is given by z and 
the phase difference by <f), it is instructive to write the geometric quantities K and r in terms 
of these physical quantities and the system parameters V,AE and A: from Eq. (5.4), 

K = 2(V 2 + (AE + A^) 2 (l - z 2 ) - V 2 (l - z 2 ) cos 2 

+2V(AE + Az)zVl - z 2 cos <p) l/2 . (5.10) 
After a short calculation K can be written as, 

hr 2 

K = 2[V 2 + (AE + Az) 2 - {H d + — ) 2 ] 1/2 , (5.11) 

where H c \ is the effective classical Hamiltonian given in Eq. (2.10), which leads to the 
integrable dynamics of z and 0. Next from Eq. (5.5) we obtain r: 

t = z(AE + Az + —= cos 4>) 
VI — z 2 
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d , (AE + Az)y/1- z 2 + Vzcos(/) ,„ n . 

+ — tan . . 5.12 

at 1/ sm 



Using the expression for H c i once again, we get 



Az 2 V 

T = H r i H 1 ; COS ( 



2 y/T 



z 



2 



at V sin 

We consider some special cases: 

(1) Interacting Bose system with no external potential: ( V = 0, A ^ 0.) 

From Eq. (2.8), setting V = 0, we get 2=constant. This in turn yields r = H c i + 
Az 2 /2=constant and K = 2(AE + Az)y/l — z 2 =constant. i.e., the underlying geometry 
is that of a circular helix with a constant pitch. 

(2) The Ideal Bose Gas in an external potential: (A = 0, V ^ 0.) 

If one considers a non-interacting Bose system then A = (see (2.9). In this limit the ki- 
netic energy term of the Hamiltonian H c \ (2.10) vanishes and hence the non-rigid pendulum 
analogy is not valid anymore. However the tunneling Hamiltonian coincides with that of a 
two-component BEC in the rotation frame approximation Here, 



H d = -V VI - z 2 cos <p + AEz (5.14) 

From Eq. (5.13), we see that, 



V ± d ^AEVT^+Vzcos^ . . 

T = H d + cos0+— tan —— — . (5.15) 

y/l — z 2 dt V sm 

Further, from Eq. (5.11), 

K = 2(V 2 + (AE) 2 -H d )^ 2 . (5.16) 

Since H c i is a constant under time evolution, Eq. (5.16) shows that the curvature K is a 
constant. However, the torsion r is time-dependent in this case. Since K is a constant, the 
path length on the unit sphere as given by the Fubini-Study metric is linearly dependent on 
time for this case. (If V and AE are made time dependent, then K is not a constant any 
more.) 

(3) The linear limit: 

For a symmetric trap with AE = in the small oscillations limit, linearizing Eqs. (2.8) in 
both z and <fi, for \z\ « 1, \<p\ << 1, we get, 

ft -< A + V >> < 5 - 17 > 
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and 

H=(A + V)^ + V^- 

This is just the harmonic oscillator limit and analytical solutions are known. The corre- 
sponding expressions for K and r can be calculated using Eqs (5.11) and (5.12). 

(4) The pendulum limit: 
For a symmetric trap with AE = 0, linearizing Eqs. (2.8) in z only, with A >> 1 we get 
the equations of a pendulum, 

dz T , . , dd . „„. 

— = -\/sin0 -r=Az 5.18 
dt Y dt y ' 

As is well known the solutions for z can be written in terms of Elliptic functions thus K and 
t can be obtained from Eq.s (5.11) and (5.12). 

Finally, for a symmetric trap AE = 0, with no linearizing approximations, though the 
analytical solution of Eqs. (2.8) can be found, it is easier to work with numerical solutions 
instead, using which K and r can be computed numerically using the expressions (5.11) and 
(5.12). 



VI. SUMMARY OF RESULTS AND DISCUSSION 

The geometric phase associated with the time evolution of the wave function of a Bose- 
Einstein condensate in a double well trap has been found using a quantum approach as in 
Sec. III. We have explicitly computed the geometric phase 4> g for both cyclic and noncyclic 
evolutions of the condensate population density difference z and phase difference <p in the 
two wells, by taking an example. In Sec IV, we have shown that the geometric phase can 
also be derived using a classical differential geometric approach, by essentially mapping the 
evolution of the two states to a framed space curve with natural moving frames along the 
curve. The unit tangent vector T to the curve has a = cos -1 z as the polar angle and 
as the azimuthal angle. As we have shown, here the geometric phase arises due to the 
path-dependent rotation of the frame perpendicular to T as the system evolves in time. 

In an experimental set up, suppose one designs a double-well trap by creating a barrier 
within a trapped condensate with N atoms, using a laser sheet. At time t = 0, let A^i(0) 
and A^O) represent the number of condensate atoms in the two neighbouring traps thus 
created, so that Ni + N 2 = N. Let the difference in the condensate densities be z(0) = 
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(iVi(O) — A^(0))/iV. Let 0(0) be the initial phase difference between the two condensates. 
We propose that in an actual experiment, immediately after creating the laser sheet barrier, 
if the density difference and the phase difference between the condensates in the two traps 
can be measured as a function of time, then by substituting these experimentally measured 
functions in Eq. (3.17), the associated geometric phase 4> g can be determined. As we have 
seen, <p g will depend on system parameters as well as initial conditions. 

Theoretically, the evolution equations for z(t) and <j)(t) are given in Eqs. (2.8a) and Eq. 
(2.8b) respectively. The trap parameters are given in Eq. (2.9). As an illustrative example, 
we have chosen the parameters 4 = 5 and AE = 0. This corresponds to an interacting 
bose system in a symmetric trap. With appropriate initial conditions (z(O),0(O)) for both 
a librational orbit and and a rotational orbit, we have solved Eqs. (2.8a) and (2.8b) for z(t) 
and <p(t) numerically. We then calculate the geometric phase <p g for the above two types of 
orbits, by using Eq. (3.17). These are plotted in Figs. 8 and 9, respectively. As is obvious, 
the last point for <p g on each of the plots, i.e., for the maximum value t = T m , corresponds 
to the geometric phase for a cyclic evolution, when the population difference and phase- 
difference between the two condensates evolve in such a way as to return to their initial 
values after a time period T m . All the other intermediate points correspond to non-cyclic 
evolutions. 

As should be clear, <p g for other parameters A and AE can also be computed by following 
our method, case by case. This would enable one to study the variation of <p g with trap 
parameters, which would be useful in designing appropriate experiments to measure it. 

The possibility of another type of experiment to study tunneling between condensates 
has been proposed by Williams et al This hinges on the fact that it has become possible 
to confine a two-component bose condensate in the same trap, as follows. Hall et al 
first trapped and cooled 87 Rb atoms in a magnetic trap in the |/ = 1, rrif = — 1 > hyperfine 
state. After condensation, it is possible to populate the \ f = 2, m/ = l> hyperfine state 
through a two-photon transition. In the presence of a weak magnetic field, these states are 
separated in energy by uo D (say). Thus two different hyperfine states can exist in the trap. 
A weak two-photon driving pulse is applied which couples the two states and consequently, 
atoms can get transferred (or "tunnel") between the two condensates. In this model, it has 
been shown |l6j| that in the mean field approximation, one obtains coupled equations for z(t) 
and <p(t) almost identical in form to Eqs. (2.8), but with A = (i.e., non-interacting) with 
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the other parameters appropriately defined for the model, and hence all our results for the 
geometric phase are applicable here as well. 

Recently, Fuentes-Guirdi et al jl^] have proposed a method for generating a geometric 
phase in a coupled two-mode Bose Einstein condensate, starting with a Hamiltonian for two 
condensates existing in different hyperfine states. In addition to the experiments of Hall 
mentioned above, condensates of 87 Rb atoms in hyperfine states |/ = l,m/ = l> and 
|/ = 2, rn/ = 2 > have been produced experimentally Q]. Likewise, condensates of 23 Na 
atoms with \ f = l,rrif = 1 > and \ f = l,m,f = > have also been created [19] . Using 
the Schwinger angular momentum representation, the Hamiltonian describing two coupled 
hyperfine states \A > and \B > can be expressed in the form \v\ 

H h f = a J z + PoJz + lo[Jx cos <f) D + J y sin <p D ], (6.1) 

Here, (J x , J y , J z ) are the components of an effective 'mesoscopic' spin J, since it can be 
shown that J is proportional to the total number of atoms N in the condensate, which is of 
the order of 10 4 . In Eq. (6.1), <pD = D t, D being the detuning frequency of the laser which 
couples the two hyperfine states. Further, 

a = {uJa - ujb) + (2J - 1)(U A - I7 fl ), A, = Ua + U ^~ Uab ^ (6 . 2 ) 

and 70 is the strength of the laser-induced drive term that couples the two levels. 
Interestingly, if we write the components of J in the form 

J = (J x , J y , J z ) = J(sina cos 0, sin a sin 0, cos a) 

in Eq. (6.1), then on setting <$£> = 0, Hhf/ J becomes identical to our Hamiltonian in Eq. 
(2.10), on identifying «o = AE, (3q = A/2 and 70 = —V. Conversely, if an external driving 
field phase (ftp is subtracted from <p in Eq. (2.10), we would essentially obtain Eq. (6.1). 
Thus our results for the geometric phase will be valid for that case too, with the appropriate 
parameters substituted. 

However, we point out that in the setting of |l^ |. one varies the BJJ system parameters 
«o and Po adiabatically to produce a closed circuit in parameter space, i.e., one considers 
an adiabatic, cyclic evolution in parameter space, with an associated geometric phase. In 
contrast, in our setting, we have fixed system parameters, for which we obtain the geometric 
phase for both cyclic and noncyclic (and nonadiabatic) evolutions on the unit sphere, i.e., 
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in "spin" space. In other words, as the population difference and phase difference between 
the two hyperfine states evolve in time, the system follows a path in this projective space, 
and as we have shown, an associated geometric phase can be defined. 

Experimental techniques to produce two condensates in close proximity has been sug- 



gested recently by Chikkatur et al 



2l( in the context of a technique to produce a continuous 



source of a Bose-Einstein condensate. It would be interesting to study the tunneling be- 
tween the condensates in such a set up if feasible, i.e., investigate the evolution of population 
differences and phase differences, to find the associated geometric phase. 

Geometric phases have been recently shown to have relevance in the implementation of 



ault-tolerant quantum computation 



221, and in the creation of vortices in a condensate 



23[. We hope that our results will have applications in these contexts as well. 
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FIG. 1: Phase portrait of BJJ evolution (Eq. (2.13)) for an interacting bose system with A = 0.5, 
in a symmetric trap. 
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FIG. 3: Phase portrait of BJJ evolution (Eq. (2.13)) for an interacting bose system with A = 5, 
in a symmetric trap. 
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FIG. 4: Phase portrait of BJJ evolution (Eq. (2.14)) for a non-interacting bose system in a 
symmetric trap with AE = 
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FIG. 7: BJJ evolution of the unit vector r (see Eq. (4.2)) on the unit sphere: Paths corresponding 
to a librational orbit and a rotational orbit (labeled r and 1 respectively in the plot) in the phase 
space portrait of the BJJ Hamiltonian for a symmetric trap with A = 5 (see Fig. 2)are shown. 
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FIG. 8: Evolution of the geometric phase as a function of time over a period for a librational orbit 
(oscillation about the zero-state) at A = 5, with initial conditions (z,4>) = (.3,0) corresponding to 
orbit 1 in Fig. 7 
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FIG. 9: Evolution of the geometric phase as a function of time over a period for a rotational orbit 
at A = 5 with initial conditions (z, <p) = (.9, 0) corresponding to orbit r in Fig. 7. 
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